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Abstract 

We study supersymmetric dark matter in the general flavor diagonal MSSM by 
means of an extensive random scan of its parameter space. We find that, in contrast 
with the standard mSUGRA lore, the large majority of viable models features either 
a higgsino or a wino-like lightest neutralino, and yields a relic abundance well below 
the WMAP bound. Among the models with neutralino relic density within the 
WMAP range, higgsino-like neutralinos are still dominant, though a sizeable fraction 
of binos is also present. In this latter case, relic density suppression mechanisms 
are shown to be essential in order to obtain the correct neutralino abundance. We 
then carry out a statistical analysis and a general discussion of neutralino dark 
matter direct detection and of indirect neutralino detection at neutrino telescopes 
and at antimatter search experiments. We point out that current data exclude only 
a marginal portion of the viable parameter space, and that models whose thermal 
relic abundance lies in the WMAP range will be significantly probed only at future 
direct detection experiments. Finally, we emphasize the importance of relic density 
enhancement mechanisms for indirect detection perspectives, in particular at future 
antimatter search experiments. 



PACS numbers: 12.60.Jv, 14.80.Ly, 95.35.+d 
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1 Introduction 



One of the hints of new physics beyond the standard model of elementary particles (SM) 
comes from the observation that most of the matter contained in the Universe is non- 
baryonic. It is therefore a desirable feature of any theory whose purpose is to extend the 
SM to provide an explanation of this unknown form of matter, commonly dubbed as Dark 
Matter (DM). 

The minimal supersymmetric extension of the standard model (MSSM) assumes the 
conservation of a discrete symmetry, called R parity PP, in order to prevent baryon and 
lepton number violating interactions, which would lead, for instance, to fast proton decay 
j2l 0]. The conservation of R parity has the nice feature of rendering the lightest su- 
persymmetric particle (LSP) stable; provided the interactions of the LSP are sufficiently 
weak, this yields an ideal candidate for dark matter. In particular, it has been shown 
that the LSP must be electrically and color neutral jl]. Among the large plethora of yet 
unobserved supersymmetric particles of the MSSM, there is a unique candidate for dark 
matter, the lightest neutralino. 

The recent analysis of the WMAP results on the cosmic microwave background tem- 
perature anisotropics, combined with other observational astrophysical data, offered a 
stringent prediction on the cold dark matter abundance (J^cdm) within the standard 
ACDM cosmo logical model 0^. The greatly improved knowledge of ^cdm, supplemented 
with accurate and complementary numerical packages for relic density and detection rates 
computations in the MSSM jUHEI, motivate a further assessment of what is known about 
neutralino dark matter in the MSSM. 

The generic structure of the soft-breaking lagrangian, which dictates most of the phe- 
nomenological features of the MSSM, is, however, still largely unknown Most studies 
devoted to supersymmetric dark matter assume either some low-energy relation between 
the values of the soft-breaking masses ^SJ^JE]; oi some underlying high-energy prin- 
ciple which organizes the soft terms at the grand unification scale or above P^ilTj. The 
main motivation for this kind of assumptions is that the number of free parameters in the 
general MSSM is huge [12,, and it is practically impossible to draw quantitative predic- 
tions out of such an unmanageable parameter space. 

The purpose of the present study is to make statistical statements about the im- 
plications of the stringent WMAP bounds on the mass and composition of the lightest 
neutralino and on dark matter detection perspectives, in the most general realization of 
the MSSM compatible with all phenomenological constraints. No a priori relations will 
be assumed among the soft breaking masses and parameters appearing in the lagrangian 
of the theory. 

The main tool we make use of is a very large random scan of the MSSM parameter 
space, from which we extract a sufficiently wide set of viable models, i.e. models which 
are consistent not only with all phenomenological constraints, but which also produce a 
neutralino relic abundance within the upper bound of the WMAP CDM allowed range. 

Though the relevant number of parameters we use in our scan is rather large (20 
overall, see Sec.|2I), supersymmetric models which feature the correct relic density within 
the WMAP range mainly fall into three categories (Sec. 

^See |Sl[7||Hlini for recent studies on supersymmetric dark matter in the light of WMAP results within 
particular SUSY scenarios. 



- higgsino-like LSP with a mass around 1 TeV. 

- wino-like LSP with a mass around 1.6 TeV. 

- bino-like LSP, with a wide mass range (up to the TeV region). 

Since in the case of a bino-hke hghtest neutrahno coannihilations or resonances are, as a 
matter of fact, unavoidable, we consider in great detail the various coannihilating partners 
and their relative effectiveness, as well as the general features of the resonant heavy Higgs 
boson s-channel exchanges (Sec. |3I). 

Finally, we study the dark matter detection perspectives in a number of search chan- 
nels: spin dependent and spin-independent direct detection, indirect detection at neutrino 
telescopes, and at antiprotons and positrons search experiments (Sec. E)). The aim of our 
analysis is to transparently compare different search strategies. To this extent, we make 
use of self consistent halo models, motivated by available numerical and observational 
data, and we present our results in the form of visibility ratios, i.e. signal-to-sensitivity 
ratios. 

We concentrate our analysis on two cases: models in which the neutralino thermal 
relic density is also consistent with the lower WMAP bound, and models which only 
fulfill the upper bound, and which, in general, have an excessively low relic density in 
the standard cosmological scenario. In this latter case, we assume the existence of relic 
density enhancement mechanisms which could affect the standard thermal relic abundance 
computations, and render models with large annihilation rates (and therefore low relic 
densities) compatible even with the lower WMAP bound. 

In this respect, a wealth of scenarios have been proposed, such as non-thermal produc- 
tion of neutralinos [TH] , cosmological enhancements due to quintessential effects ,20] , 
to anisotropic cosmologies producing an effective shear energy density |221; or to 
scalar-tensor theories [231 • Under the assumption of the existence of such enhancement 
mechanisms, no rescaling procedure is in order, for models underproducing dark matter 
in the standard cosmological scenario of thermal production. Indirect detection rates, 
which essentially depend on the same annihilation rates determining the neutralino relic 
abundance, will then be considerably larger than for model whose relic density falls within 
the WMAP range. 

We point out that current data on dark matter detection do not significantly constrain 
the viable supersymmetric parameter space, at least in a statistical sense. Regarding fu- 
ture perspectives, spin-independent direct detection is the most promising search strategy 
for models whose neutralino relic abundance lies in the WMAP range. We also find that 
neutrino telescopes have a limited parameter space reach, and can only probe models 
with relatively low masses. We however point out the correlations and complementar- 
ity between neutrino telescopes searches and direct dark matter detection. Finally, we 
emphasize that models with large annihilation rates, assuming some relic density enhance- 
ment mechanism, give spectacular signals at antimatter searches, and are even already 
constrained by available data. 



2 MSSM Parameter Space Scan 



2.1 Motivation 

Supersymmetry, if it exists, must be spontaneously broken. The precise mechanism of 
supersymmetry breaking is, however, not known. From a practical point of view this 
difficulty is overcome by introducing into the MSSM effective Lagrangian extra soft terms 
which break supersymmetry explicitly. The collection of such terms is known as the 
soft-breaking Lagrangian (see ^2] for a recent review), and reads 

-^rh\^h\hi + m^2^2^2 + {Bfihih2 + h.c.) 

+ {A^Y^h,d^qL, + A^^h^,, + A^Y^^h^L, 

+ ^M,BlBl + ^M^WIWI + ImsG^G"^ + h.c). (1) 

which respectively corresponds, line by line, to squark masses, slepton masses, Higgs 
masses, trilinear couplings, and gaugino masses. The most general Csoft introduces a 
huge number of free parameters (more than 100), which makes any phenomenological 
study based on a blind approach to the MSSM parameter space very problematic. The 
usual way of studying supersymmetric effects is to assume specific frameworks in which 
all soft-parameters depend on a few inputs given at a certain high energy scale. Among 
them, minimal supergravity (mSUGRA) models |^ have received particular attention. 
mSUGRA models are defined in terms of four continuous parameters and one sign: 

tan/3, mo, Oq, mi/2, sign(/i) , (2) 

which determine the whole set of soft breaking terms at the unification scale Mqut through 
the relations: 

Mi = mi/2, {m'j)ij = ml5ij, ml^ = fhl^ = ml, Af = oq. (3) 

Thus, at Mgut all scalars have a common mass mo, and all gauginos feature the same 
soft mass mi/2. 

The allowed parameter space of mSUGRA models, compatible with all phenomeno- 
logical and cosmological constraints have been determined in a number of works (among 
recent studies see e.g. jHl EH])- It turns out that the most important constraint comes 
from the requirement that the neutralino relic density does not exceed the relic density 
of cold dark matter. 

In mSUGRA, the lightest neutralino is usually a bino-like neutralino, which has a 
small annihilation cross section, and therefore tends to produce a large relic abundance. 
Indeed, apart from a small region at low neutralino mass, only three regions fulfill the 
WMAP upper bound on the relic density: 

• The stau coannihilation strip: Along this region the stau is almost degenerate with 
the lightest neutralino, and coannihilation processes help suppress the relic density. 



• The Funnel region: In this case there is a resonant enhancement of the bino-bino 
annihilation cross section through s-channel heavy Higgs bosons exchange. A nec- 
essary condition required for such a resonant enhancement is a large value of tan /5, 
needed to fulfill the relation 2m^ ^ m^. 

• The focus point region: These are narrow regions with very large values of mo 
yielding, through electroweak symmetry breaking conditions, a low value of the n 
parameter, and implying a lightest neutralino with a non-negligible, or even a large 
higgsino fraction. 

A legitimate question one can ask at this point is, therefore, whether or not mSUGRA 
is a benchmark scenario for SUSY dark matter with respect to the bulk of the general 
MSSM parameter space, or if it is a somehow theoretically biased framework. In this paper 
we try to give an answer to this issue too, taking a different approach to the study of 
supersymmetric models. Instead of choosing a particular model of SUSY breaking, we 
carry out a statistical analysis of a huge number of models with randomly generated soft 
parameters. We do not make any simplifying hypothesis for the low energy parameters, 
and we therefore do not assume any scalar universality or gaugino unification relations. 
In this way we wish to get a feeling of the typical predictions of SUSY models, and to 
compare them against mSUGRA or any other specific scenario. 



2.2 The Scan 

Many of the parameters of Csoft are severely constrained because they would imply FCNC 
or CP violating effects at a rate which is already ruled out by the experiments. For 
simplicity, we will assume, as usually done in the literature, that all soft parameters are 
real, so that supersymmetry breaking does not introduce new sources of CP violation. 
To suppress potentially dangerous FCNC, we will set to zero all off-diagonal elements 
in the sfermions masses, and assume that the first and second generation of squarks are 
degenerate (for a recent discussion see ^3). Squarks masses therefore have the form 




with rriq, mg^, and m^g arbitrary numbers. Slepton masses, on the other hand, 
contain three independent entries each. 

The trilinear couplings At, A^, Aj. and A^ are allowed to have both signs. All other 
trilinear couplings are neglected^. Gaugino masses are independent of one another, and 
negative values for M2 are also considered. The remaining low energy parameters we take 
into account are tan/?, /i and m^, the mass of the pseudo-scalar Higgs boson. 

We scan the resulting 20-dimensional parameter space using a uniform probability dis- 
tribution, tan /? takes values between 2 and 50, whereas all mass parameters are generated 
in the interval (50 GeV, 5 TeV) possibly with both signs, according to Table ^ From this 
set of low energy parameters, one can determine the mass spectra and mixing matrices of 
the superparticles. 

^We include for its relevance in the computation of (g — 2)^. 



Name 


# of parameters 


Symbol 


Range 


Gluino and Bino Masses 


2 


Ml, Ms 


(50 GeV, 5 TeV) 


Wino Mass 


1 


M2 


± (50 GeV, 5 TeV) 


Left-handed slepton masses 


3 


mr 

I'i 


(50 GeV, 5 TeV) 


Right-handed slepton masses 


3 




(50 GeV, 5 TeV) 


1st and 2nd family squark masses 


1 


ruq 


(50 GeV, 5 TeV) 


3rd family squark masses 


3 




(50 GeV, 5 TeV) 


Third family trilinear couplings 


3 


-^ti -^hi 


(-5 TeV, 5 TeV) 


Pseudo-scalar Higgs boson mass 


1 




(50 GeV, 5 TeV) 


Muon trilinear coupling 


1 


A. 


(-5 TeV, 5 TeV) 


/i 


1 




±(50 GeV, 5 TeV) 


tan (3 


1 


tan P 


(2, 50) 


Total 


20 







Table 1: List of the MSSM parameters taken into account in our scan, and the range on 
which they are allowed to vary. 

All the results we present in this paper are obtained from this kind of procedure. Had 
we chosen to scan over a smaller number of parameters (assuming additional relations), or 
to drastically change the range on which they vary, the results might have been different. 
We claim however that this scan covers what can be considered a natural parameter space 
range, and that it is largely free of theoretical prejudices. 

2.3 Phenomenological Constraints 

We apply the following phenomenological constraints: 

• The spectrum: The presence of non-zero trilinear couplings could give rise to tachy- 
onic sfermions, so the first consistency requirement we ask is to exclude all such 
unphysical models. 

• Neutmlino LSP: In the MSSM the lightest supersymmetric particle is stable and 
therefore must be an electrically neutral and not strongly interacting particle jl]. 
Since the sneutrino, in the MSSM, has been shown not to be a suitable dark matter 
candidate j^, the only possibility we are left with is the lightest neutralino. We 
therefore also require the LSP to be the lightest neutralino. 

• h ^ 57: In SUSY models with minimal flavor violation the decay b s'-f proceeds 
through the tW and tH~^ loops, in addition to the SM contribution from the tW 
loop. The branching fraction BF{b 57) has been measured by the BELLE, 
ALEPH, and CLEO collaborations. A weighted averaging of these measurements 
of -B — ^ Xs7 decays at CLEO and BELLE lead to bounds on the branching ratio 
b ^ s'-f. We will require 

2 X 10^^ < BR{b 57) < 4.6 x 10"^ (5) 

The calculation of BR{b 57) is carried out with the latest release of the micrOMEGAs 
package j2Hl; which includes an improved NLO and the charged Higgs contributions 
as well as a beyond leading-order treatment of large tan (3 effects. 



Direct Accelerator Searches: We also take into account the limits derived from the 
unsuccessful searches for supersymmetric particles and the Higgs boson. In Table |21 
we detail on the mass limits we impose on each SUSY particle. 



Name 



Mass Limit (GeV) 



Charginos 
sneutrinos 
charged sleptons 
sbottoms 
stops 

other squarks 
gluino 
Higgs boson 

Pseudoscalar Higgs boson 



103.5 
43.0 
95.0 
91.0 
86.4 
100.0 
195.0 
114.0 
85.0 



Table 2: Mass limits 



• Relic density: Finally, we require that the neutralino relic density lies below the 
WMAP 2 - cr upper bound jS] , 



The WMAP lower bound (I^dm^^ > 0.09) is not imposed since it may well be that 
neutralinos are only a sub dominant dark matter component. For the evaluation of 
the neutralino relic density we again make use of the program micrOMEGAs, which 
takes into account all tree level annihilation processes as well as all possible coan- 
nihilation processes in the MSSM. We also cross-checked our statistical results for 
a sample of models with the DarkSUSY package, and verified that they are overall 
consistent with each other. 

Although our main focus is on supersymmetric Dark Matter, we will first use our scan 
to quantify the effectiveness of the different constraints mentioned above. The results are 
illustrated in Fig. ^ Different colors indicate different samples: Red for the whole sample. 
Blue for those models surviving the red constraint, and green for those surviving the blue 
one. We see that 10% of the models contain a tachyon; of the remaining models, 20% 
are ruled out by the 6 — > S7 constraint, 25% are excluded by direct searches and almost 
80% contain a LSP different from a neutralino. Finally, 30% of the models which are not 
excluded by any of the previous tests give rise to a neutralino relic abundance larger than 
the WMAP upper bound. The surviving models will be called viable, and are the only 
ones we will consider from now on. Our sample consists of 10^ viable models. 

2.4 A statistical analysis of 6 ^ 57, nih and {g — 2)^ 

We take at this point the opportunity to assess the statistical results concerning three 
quantities sensitive to supersymmetric effects, i.e. the rare decay b — >■ sj, the mass of the 
lightest CP-even Higgs boson and the muon anomalous magnetic moment {g — 2)^. 

In panel (a) of Fig. |21 we show a histogram with the value of BF{b 57) in units 
of 10""^. Notice that most models predict a value close to the central value predicted by 



QoMh^ < 0.13. 



(6) 



Tachyons b->SY Searches LSP q ^2 



Figure 1: Effectiveness of the different constraints imposed on supersymmetric models. 
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Figure 2: (a): BF{h 57) in units o/lO"'*. (b): 
bound mil ^ 132 emerges from our scan. 
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the SM 3.4). That is, the supersymmetric contributions to 6 — > 37 tend to be very 
small, typically signaling a heavy SUSY spectrum. On average, the SUSY contributions 
are found to be positive, and to increase, therefore, the SM value of BF{h — > S7). 

Among the different limits shown in Tab. El the Higgs Boson mass plays a special 
role. Indeed, one of the generic predictions of the MSSM is a light, "standard-model-like" 
Higgs boson. The computation of the Higgs mass is carried out with the FeynHiggsFast 
package, which includes radiative corrections up to the two-loop level. For the top mass 
we set uit = 175 GeV. In Fig. |21(b) we show a histogram of the Higgs boson mass in units 
of GeV. Notice that there is a smooth peak around rrih ~ 119 GeV, and that no models 
were found with ruh > 132 GeV (this limit is very similar to what reported in earlier 
analysis, see e.g. [22] )• 

Due to a long standing uncertain theoretical and experimental situation, we did not 
take into account any bound deriving from the anomalous magnetic moment of the muon. 
Nevertheless, we show in Fig. |3] a histogram of the supersymmetric contribution to = 
{g — 2)^/2. The calculation is done with the routine provided with micrOMEGAs. It is 
clear from the figure that most models predict a small value of (5a^: 

-2.0 X 10"^° < 6a^ < 2.0 x 10"^° . (7) 
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Figure 3: The supersymmetric contribution to the anomalous magnetic moment of the muon 
(a^ = {g- 2)^/2) in units of 10'^° . 

Using the latest experimental and theoretical results, based on data from e'^e~ — > hadrons, 
the deviation Sa^^ was found to be pH] : 

Sa^ = (27.1 ± 9.4) x 10"^° (Davier et al.) [30j (8) 
5a^ = (31.7 ± 9.5) X 10"^° (Hagiwara a/.) ^ (9) 

If, instead, the r decay data is used to determine the hadronic vacuum polarization, the 
deviation is smaller: 

5a^ = (12.4 ± 8.3) x 10"^° (Davier et al) [SD! ■ (10) 

In a conservative approach, such as that outlined in Ref. (22], we can therefore conclude 
that the bounds stemming from a^ only constrain a marginal portion of the viable MSSM 
parameter space. 

3 The Nature of Neutralino DM after WMAP 
3.1 Relic Density of Neutralino DM 

The relic density of the lightest neutralino depends critically on two factors: the neutralino 
mass rriy and the neutralino interactions. In the MSSM, the lightest neutralino is a linear 
combination of the gauge eigenstates, which are dubbed bino, wino and higgsino; its 
interactions are therefore determined by the relative bino-, wino- and higgsino-content. 
We show in Fig. |3] a plot of the relic density as a function of the neutralino mass for 
a bino-, a higgsino- and a wino-like neutralino. All other superparticles are assumed to 
have a mass much larger than and given by msusY ~ ct ■ (Notice that doing so, 
we are neglecting the possibility of coannihilations or resonances through an s-channeF). 
From the plot we see that the dependence with a is particularly relevant for a bino-like 
neutralino. Some general conclusions can be drawn from this plot: 

■^Assuming that the scalar masses he well above the gaugino masses amounts to resorting to a scenario 
similar to that of the finely tuned MSSM j^Sl, sometimes also dubbed Split Super symmetry, see |34) . 



• A bino-like neutralino typically produces a relic abundance that is much larger than 
the WMAP bound. In fact, we will see that viable models with bino-like neutralinos 
require the presence of peculiar mechanisms which suppress the relic density. 

• A higgsino-like neutralino must have a mass of the order of 1 TeV in order to produce 
a relic abundance consistent with the preferred WMAP range. 

• A wino-like neutralino generates a relic abundance in the WMAP range only pro- 
vided ~ 2 TeV. 

In Fig. |S1 we show a histogram of the total neutralino relic density, and of the relic 
density of bino-, wino- and Higgsino-like neutralinos. A neutralino is said to be bino-like if 
its bino component is larger than its wino and higgsino components. Wino- and higgsino- 
like neutralinos are defined analogously. From the first diagram, we learn that most models 
produce a relic density that is much smaller than the WMAP bound (the gray vertical 
band), with a moderate peak around Q^h"^ ~ 10^^. This picture is clearly very different 
from the usual situation in supergravity inspired frameworks, where the neutralino relic 
abundance is typically found to be excessively large. The second diagram reveals that 
bino-like neutralinos typically yield a large relic abundance, and that they contribute 
significantly to the WMAP range. Models with wino- and higgsino-like neutralinos, on 
the other hand, feature a relic density well below the WMAP range, as apparent from the 
Wino and Higgsino histograms. 

3.2 Neutralino Mass and Composition 

The composition of a typical neutralino markedly depends on whether or not we demand 
the WMAP lower bound to be fulfilled. This fact is illustrated in Fig.lHlwith pie diagrams 
showing the relative number of bino-, wino-, and higgsino-like neutralinos. In (a), where 
the lower bound was not imposed, the lightest neutralino is predominantly a higgsino or 
a wino, with only a small fraction of binos. In (6), where the relic density is required to 
lie within the WMAP range, the lightest neutralino is mostly a higgsino, with a sizeable 
fraction of binos and a moderate contribution from winos. 

In Fig. 13 we show histograms of the lightest neutralino mass for all viable models. 
As expected from statistical reasons, the probability of having a neutralino LSP with a 
given mass m decreases with m. Notice that the mass probability distribution P{m~^) of 
a wino (and of a higgsino) follows the statistically expected functional form, dictated by 
the requirement that M2 (respectively /i) is the lightest among np ~ 15 mass parameters 
which range in the [^up^ ""^down] interval, i.e. the following power-law: 



On the other hand, the distribution we find for the bino masses reflects the one obtained 
under the requirement of having at least one coannihilating partner, whose mass lies 
between m-^ and (1 + S)m^, i.e. 



P{mx) oc (mup - m^) 



(np-l) 



(11) 



P (m^) oc (mup - m^) 



(np-l) 
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Figure 4: The neutralino relic density as a function of for a bino-, a wino- and a higgsino- 
like neutralino. All other superparticles are assumed to have a mass much larger than m^: 
mrisusy ~ o • ""^x' '"'^^^ a = 3, 5, 10. 

which, in the hmit of small 5, features a maximum at 



in reasonable agreement with what we statistically find in the case of binos^. This is 
clearly a first strong indication of the (statistical) relevance of coannihilations for bino- 
like neutralinos; although coannihilations may appear as a finely-tuned mechanism of relic 
density suppression, cosmology seems to point, instead, to their "naturalness". We will 
carry out a more detailed analysis of this point in Sec. 0] 

In Fig. ISl we show histograms of the neutralino mass for models with a relic density 
in the WMAP range. We see that most models have a neutralino mass in the range 
800 GeV< rriy^ < 1200 GeV. Binos are usually light (m^ < 1 TeV) whereas winos and 
higgsino masses cluster around ^1.6 TeV and 1 TeV. Since the LHC reach for 

higgsino-like or wino-like neutralinos will likely be well below an LSP mass of 1 TeV, we 
point out that SUSY models with a relic abundance lying within the WMAP range and 
featuring a wino or higgsino-like lightest neutralino, will not be detectable at the LHC. 

4 Relic Density Suppression Mechanisms 

In the previous sections we pointed out the well known fact that the neutralino relic 
density critically depends on the composition of the LSP in terms of its bino, wino and 
higgsino components. The annihilation cross section of a bino is by far smaller than those 
of a wino or a higgsino: first, the couplings involving a bino are smaller than those, for 

*The location of the maximum in the histogram is at sUghtly lower masses due to the presence of 
non-coannihilating light binos. 



= mup/{np — 1) ^ 350 GeV, 



(13) 




Figure 5: Histograms of the neutmlino relic density. The first panel shows the statistical distri- 
bution of the neutralino relic density for all models. In the other panels, the relic density of a 
bino-, a wino- and a higgsino-like neutralino are shown. The gray vertical band corresponds to 
the WMAP range. 




Figure 6: Pie diagrams showing the composition of the lightest neutralino. In (a) we include 
all viable models; in (b), only those with a relic density within the WMAP range. 

instance, of a wino, since gi < g2] furthermore, the number of final states in which a 
bino can annihilate is much smaller than that of winos and higgsinos, and the latter can 
directly annihilate, without the need of an intermediate supersymmetric partner, into, for 
example, a couple of gauge bosons, whereas binos cannot. This translates into a critical 
dependence of bino annihilation cross section on the supersymmetric spectrum, which is. 
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Figure 7: Histograms of the neutralino mass. The first panel shows the mass of the neutralinos 
for all viable models. In the others, we examined each kind of neutralino separately. 

on the other hand, absent in the case of winos and higgsinos. Last but not least, the 
mass matrix structure of charginos and neutrahnos imphes that a wino-hke LSP has a 
quasi degenerate hghtest chargino, yielding a large coannihilation contribution; for the 
same reason, higgsinos are quasi degenerate both with the lightest chargino and with the 
next-to-lightest neutralino, again with large corresponding coannihilation effects. In this 
respect, the features of Fig. HI come not as a surprise: the relic density of higgsinos and 
winos has a weak dependence on the SUSY particle spectrum, and is orders of magnitude 
suppressed with respect to that of binos. 

The bottom line is, on the one hand, that the mass spectrum of higgsinos and winos 
whose relic density falls within the WMAP range clusters around definite values of the 
neutralino mass, the only relevant parameter. On the other hand, we will hereafter show 
that, in most cases, binos are compatible with the dark matter abundance only if the 
particle spectrum is such that coannihilations or resonant annihilation channels are open. 

To make the previous statement more quantitative, in this section we limit our dis- 
cussion to neutralinos with a bino purity larger than 0.9. We find that: 

1. Among the binos compatible with the upperWViKP bound, 91% have at least one 
coannihilating partner^ and 15% have a resonant annihilation cross section with the 

particle is considered a coannihilating partner if the ratio Afi/il between the difference of the relic 
density computed without taking into account coannihilations and that with coannihilation f2, over is 
larger than 1% (see next section for details). 
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Figure 8: Neutralino masses giving a relic density in the WMAP range. 



heavier Higgs boson mass ^ 2 ■ m^. Only 2% of binos do not feature any relic 
density suppression mechanism. 

2. Binos which produce a relic density within the WMAP range have a coannihilating 
partner in 86% of considered cases and a resonant cross section in 18% of cases. 4% 
of binos does not feature any relic density suppression mechanism. 

3. The largest bino mass compatible with the WMAP bound, and such that neither 
coannihilations nor resonances are present, is around 160 GeV. 

It is therefore clear that in the large majority of cases, and for all masses larger than few 
hundreds GeV, binos can be dark matter candidates in the general MSSM only if either 
coannihilations or resonances are present. In this respect it is worthwhile to analyze in 
greater detail both relic density suppression mechanisms. 



4.1 Coannihilations 

The importance of co- annihilations, i.e. cases where the number density of a species 
depend not only on its own annihilation cross section but also on the effect of the annihi- 
lations of another species, was first recognized in the seminal papers of Ref. P31 l36j. When 
the mass of the stable LSP is close to the mass of other particle species, the evolution of 
the number density of LSP's, nLSP, and that of the other particles, ni=i^,,,^N, are tightly 




Figure 9: Effectiveness of bino coannihilations in the MSSM: a plot of Ail./ Q as a function of 
the mass splitting between a bino-line neutralino and all possible coannihilating partners in the 
MSSM. 

correlated, and must be tracked together while solving the Boltzmann equation 

^ = -3Hn - (aeflff ) (n^ - n^^) , 
where the number density is now given by 

Af 

n = riLSP + y^^rij. (15) 

i=l 

In the equations above H denotes the Hubble constant, (acsv) is the thermal average of 
the effective cross section times the relative velocity, and rieq is the equilibrium number 
density, which, to a good approximation, may be taken to be equal to the Maxwell- 
Boltzmann thermal distribution. The main effects of taking into account coannihilations 
are two-folds: on the one hand one modifies the effective interaction cross section; on 
the other hand one alters the number of degrees of freedom which enters in the game. 
Therefore, if the cross section of the coannihilating partner is much more efficient than 
that of the stable species, the net result will be a reduction of the number density nLSP- 
In the opposite case, i.e. when the extra coannihilating degrees of freedom carry less 
efficient annihilations, the outcome can be an increase of the final asymptotic nLSP- Since, 
however, the coannihilating partners have typically non-zero electric or color charges, 
while neutralinos are neutral, this second effect is rather unusual (for exceptions see jHI 
in the mSUGRA context and jHZI in the gluino coannihilation model). 

Since the discovery of the generic mechanism of coannihilations, many dedicated stud- 
ies have analyzed the impact of considering various coannihilating partners. In particular, 
within the framework of the constrained MSSM, the NLSP is found to be, in the low tuq 



(14) 



region, the lightest stau [SHI EHl HOI E] • In the focus point region of the CMSSM, instead, 
a non-trivial higgsino fraction may give rise to chargino and next-to-lightest neutralino 
coannihilation processes, see e.g. [HI |321 1321 El • Again within the CMSSM, at large scalar 
trilinear couplings, it may well be that the next-to-LSP is the stop, whose coannihila- 
tions were considered in |S1 1^ E] • Minimal deviations from the assumed universality 
of scalar soft-breaking masses (SBM) have been shown to lead to other viable coannihi- 
lating partners: for instance, lowering the scalar SBM of the particles belonging to the 
5 representation of SU{5) gives raise to coannihilations with the lightest bottom squark 
and with the tau sneutrino ^Hl El ■ Relaxing the assumption of universality at the high 
energy (grand unification) scale in the Higgs sector may also give rise to sneutrino and 
other coannihilating partners which are not present in the CMSSM |3HJ|311. Finally, in 
a recent analysis it has been shown that the strongest coannihilation processes in the 
MSSM are those with the gluino, and that the maximal bino mass is reached precisely in 
the gluino coannihilation tail [HTj . 

The effect of including coannihilations in the computation of the relic density of neu- 
tralinos has been carried out for particular coannihilating partners, see e.g. in Ref. 0. 
In jnZl we analyzed the relevant case of bino coannihilations with all possible partners, 
assuming, within an effective MSSM, that all relevant soft breaking masses are three times 
larger than the bino mass mi, except for the particular coannihilating partner mass, which 
was taken to be close to mi. The relevant parameter we used was the relative splitting 
between the bino mass and the mass of the coannihilating particle mp, 

nip — rriy , , 

Am = — ^ ^ 16 

m^ 

The results are mostly independent of both the absolute size of mi and of the details of 
the spectrum of the other SUSY particles. Moreover, the dependence on tan was found 
to be not critical. 

We show in Fig. El the relative difference in relic density without and with coannihila- 
tions ini 

Af2 coann ^coann 

= — . (17) 

A. 

The lower limit on Afi/f2 has been set to 1%, since this is the typical numerical accuracy 
of the numerical packages we employed to carry out the computations jHl EHl • The plot 
was taken at a neutralino mass m^ = 400 GeV and at tan/3 = 30. 

From the plot we clearly deduce that gluino coannihilations are the strongest possible 
bino coannihilation processes in the MSSM, as it might be expected considering the strong 
gluino-gluino annihilation cross section and the number of final SM states, much larger 
than that of squarks. Binos can also annihilate with a quasi degenerate wino without 
altering its bino purity (the same does not hold true for the higgsino, i.e. in case mi ~ /x, 
due to the neutralino mass matrix structure), and therefore undergo chargino and next- 
to-lightest neutralino coannihilations. These processes are found to be of the same order 
of magnitude of those involving up and down squarks (or two of them at the same time). 
Much more suppressed are, instead, slepton coannihilations. 

Fig. El allows us to make a quantitative statement about models where coannihilations 
are or not relevant in reducing the relic abundance of binos to acceptable levels. We will 
declare that a model has a coannihilating partner P if the mass splitting of that partner 
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Table 3: Mass splitting within which coannihilations with different particles are effective. 

from the bino mass, Amp, is such that the corresponding AQ/Q is larger than 1%. In 
this respect, we show in Tab. |3]the cuts on the mass splitting corresponding to all possible 
bino coannihilating partners, as well as the symbols we will use in the following figures. 

We then consider all the models in our scan featuring a bino-like LSP with a purity 
larger than 90%, and whose relic abundance falls within the WMAP bounds, and look 
for the number of models exhibiting a given coannihilation channel. Our results are 
summarized, with the same symbols as in Tab. El in Fig. (a ). We see that in most 
models with the correct neutralino relic density a squark coannihilation is present: this is 
mainly due to the fact that squark coannihilations are strong, but also to the number of 
different particles falling into this category (e.g. 6 up-type squarks). Therefore, to single 
out the relevance of coannihilation processes with a given species, we divide the number of 
models characterized by the occurrence of a given coannihilating partner by the number 
of possible coannihilating partners of that given species. We show our results in Fig. 
{b ): remarkably, we obtain the same hierarchy of coannihilation efficiency outlined in 
Fig. El gluino coannihilations are the strongest one, followed by squark coannihilations, 
chargino and NL-neutralino and finally sleptons. 
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Figure 10: Bino- coannihilations in the MSSM. In [a), we show, for each possible coannihilating 
partner, the percentage of viable models with bino-like neutralinos in which coannihilations with 
one of those partners take place. In (6) we have normalized the data from (a) dividing by the 
number of particles in each class (e.g. by six for I, by three for V, etc. ). 




Figure 11: The relic density in the presence of a resonant annihilation cross section, as a 
function of ArriA, for different values of tan P and different neutralino masses (a). In (6) we 
display an histogram of Atua for all viable models with a bino-like neutralino and in the 
WMAP range. 



4.2 Resonances 

The computation of the rehc density of a thermally produced particle x needs particular 
care when the mass of x li^s close to half the mass of another particle P in which x can 
annihilate through an s-channel P exchange. In this case, in fact, the propagator of the 
s-channel diagram has a resonant behavior. This is the case, for instance, in the MSSM 
for the heavy CP-odd A and CP-even H Higgs bosons^. The calculation of the thermally 
averaged cross section for the resonant XX annihilation cross section needs the thermally 
averaged treatment of a cross section of the Breit-Wigner form see e.g. [S^ . 

The effects of a heavy Higgs resonance on the annihilation cross section of neutralinos 
depends on three inputs: first, the mass splitting between the annihilating particle x 
and half the mass of the resonantly exchanged particle A, which we quantify through the 
parameter 

AmA = ^; (18) 

2 ■ 

second, the total decay width of the resonance Va and third its mass |52j. Finally, the 
ratio of the naively computed cross section and the correct one depends on the relevant 

^Another possibility is a resonance with the lightest CP-even Higgs, which would however require a 
neutralino mass below 60-70 GeV: this case has therefore not been considered in the present study. 



XX^ couplings and on the couplings of the resonance to the final particle states. For 
recent studies dedicated to the effects of resonant annihilation cross sections see [S3]- 

We show in Fig. ^2 the overall non-trivial combination of the effects of these input 
parameters. In the upper panel we show the dependence of the neutralino relic density 
fi^/i^ on Am^ at fixed = 200 GeV for three different values of tan/3, while in the 
lower panel we fix tan /? = 30 and analyze three different m^=200, 500 and 1000 GeV. The 
plot clearly highlights the effect of the resonance around Am^ = 0, and in particular the 
results of the thermal averaging: in fact, the maximal effects are at positive ArriA ~ 1%, 
since particles with mass < ?7t,^/2 in general have larger thermal energy contributions 
than those with > rriA/^. 

In the upper panel we see the effect of a broader resonance (the width of rriA grows 
with tan/5) on the relic density: for larger tan/3 = 50 resonance effects are effective for 
a larger Am^ range. In the lower panel, instead, it is shown how larger masses reduce 
the efficiency of resonance effects and shrink the Am^ range where resonance effects are 
relevant. 

Resonances have drastic consequences on the cosmological constraints on the MSSM 
parameter space, mainly when the XX annihilation cross section is relatively low, as it is 
the case for a bino-like LSP. The plot on the right shows the overall statistical effects of 
the A pole resonance showing the fraction of binos in a given bin of Am a- We clearly see 
that the privileged values are those around Am^ = 0, as expected. 

5 Direct and Indirect DM detection 

The present section is devoted to the discussion of direct and indirect dark matter de- 
tection signals in the general MSSM (for recent reviews on the subject of dark matter 
detection see The issue of the comparison between different search strategies, and 

of the complementarity among them, has since long been the subject of various investiga- 
tions (see e.g. [SHI)- the present context we propose an approach which features three 
major novelties: 

1. We adopt self-consistent halo- models, along the lines of Ref. where two ex- 
treme benchmark scenarios for the dark matter distribution in the Milky Way have 
been outlined, both being soundly motivated from the point of view of structure 
formation, and consistent with all observational data available and with numerical 
simulations. The two models represent two different back-reaction mechanisms of 
the baryon infall on the dark matter distribution, giving rise to a non-cuspy profile 
(the i?ur/i;ert profile) and to a steeply cuspy profile {Adiabatically contracted profile). 
The Burkert profile is a conservative halo-model of the "cored" type: the central 
cusp gets smoothed out through large angular momentum transfer between the 
baryonic and the dark matter components. In the adiabatically contracted profile, 
optimistic from the point of view of dark matter detection, the starting dark matter 
distribution (assumed to be the CDM profile of Ref. ^7\) undergoes an adiabatic 
contraction with no net angular momentum transfer: The CDM cusp therefore gets 
increased, approaching a rather steep density profile towards the Galactic center 
||56j . We stress that both profiles satisfy all dynamical constraints, and that the 
velocity distributions are computed self-consistently. Moreover, they represent two 
extreme cases: it is reasonable to assume that the actual halo profile of the Milky 



Way lies in between these two cases. In this respect, our results can be regarded 
as representative lower and upper limits on SUSY dark matter detection rates. Let 
us stress that, in order to compare different search strategies, it is mandatory to 
have consistent profiles for dark matter densities and velocities, and both must be 
compatible with observational data. 

2. We make use of a statistical analysis based on Visibility Ratios, i.e. signal-to- 
sensitivity ratios, which allows a transparent comparison between different methods 
and between current and projected sensitivities. Our aim is to estimate the statisti- 
cal relevance of a given search strategy. We claim that this approach may be more 
quantitative and helpful, in the evaluation of the perspectives of detection in the 
most general SUSY setup, than usual scatter plots, which hardly bring any quan- 
titative statistical information. We will however also make use of scatter plots in 
order to show correlations among different quantities, and to outline the dependence 
of detection rates on the neutralino mass, composition and relic abundance. 

3. In view of recent results about the occurrence of relic density enhancement mecha- 
nisms in modified cosmological contexts ^3 1201 1^ 122 123 ? or through non-thermal 
neutralino production [iH] , we will not apply any rescahng procedure in the compu- 
tation of detection rates. Models with a low relic abundance may well be responsible 
for the whole inferred amount of dark matter in the Universe, and provide, thanks 
to large neutralino annihilation amplitudes, significant indirect detection rates. We 
will nevertheless always single out, for definiteness, the case of models which provide 
the required neutralino abundance in the standard thermal cosmological scenario. 

5.1 Spin Dependent and Spin Independent Direct Detection 

We show in Fig. the results of the scan for direct detection experiments in models 
with a relic density within the WMAP range. On the left we plot the spin-independent 
neutralino-proton cross section (results for the neutralino-neutron cross sections are very 
similar) including, in black, the current limits from the Edelweiss experiment [58j and, in 
green, the projected future sensitivity of the Xenon 1-ton facility ISHl, computed for the 
Adiabatically contracted halo profile^. The code for the scatter plot is as follows: black 
circles refer to bino-like neutralinos, red squares to higgsino-like and blue diamonds to 
wino-like. 

The right part of the figure shows instead our results for the spin- dependent YieniiaXmo- 
proton cross section, with current constraints from various existing experiments indicated 
with a black line and the future projected NAIAD sensitivity in green [HT]. As evident 
from the figure, spin-dependent searches always feature an experimental sensitivity which 
is unable to probe any viable SUSY model in our scan. 

Fig. El collects the statistical analysis of our results for spin-independent searches, 
as a function of visibility ratios, for the adiabatic halo profile. The histograms to the 
left show the cumulative results for all models, while those to the right refer to models 
within the WMAP range (hence those shown in the scatter plots of Fig. [T^ . The two 
upper histograms refer to the current limits (those indicated with black lines in Fig. IT^. 

^In the case of direct detection and neutrino telescopes, the dependence on the halo profile, which is 
merely local, is rather mild, see |56| . 




Figure 12: A scatter plot of the Spin Independent (a) and Spin Dependent (h) neutralino-proton 
scattering cross sections for models whose relic density lies within the WMAP range. Black 
circles represent models in which the lightest neutralino is predominantly a Bino, red squares 
represent Higgsinos and Blue diamonds Winos. The black lines indicate the current exclusion 
limits \E3i, while green lines the projected exclusion limits at future experiments p^ . 
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Figure 13: A statistical analysis of the Visibility Ratio (signal-to-sensitivity) for the Spin In- 
dependent neutrahno-proton cross section for all models with a sufficiently low relic abundance 
(plots to the left) and for models with a relic abundance in the WMAP range (plots to the right). 
The two upper plots refer to the current exclusion limits, as in Fia. (a), black line, while the 
two lower plots to the future exclusion limits, Fig.^^(a), green line. 



and the two lower ones to future exclusion limits (green lines in Fig. [T2|) . The vertical 
blue lines indicate the visibility threshold: models lying to the right of the blue line are 
above current or future sensitivities, those to the left currently give no (or will not give, 
at future experimental facilities) detectable signals. While a negligible fraction of models 
lie above current exclusion limits, approximately 20% of models (in the WMAP range or 
with lower relic abundances) will be accessible to Xenon 1-ton jSH] , or to similar detection 
experiments. 

A few comments are in order: first, spin dependent rates are found to be at least one 
order of magnitude below future sensitivities^. The potential of spin dependent searches 
is then statistically negligible. Second, current exclusion limits from spin-independent 
searches are still far from probing a significant portion of the general MSSM parameter 
space; on the other hand, future experiments will be able to probe 20% of viable models. 
The scatter plots also show the mass clustering of pure higgsinos and of winos around, 
respectively, 1 TeV and 1.6 TeV: the dependence of the neutralino-proton scattering cross 
section on the details of the SUSY spectrum, particularly concerning the Higgs sector 
(the sign of the /i parameter and the masses of the CP-even h and H neutral Higgses), 
and, to a less extent, the squark sector yield a scatter in cr^p which can well be as large 
as four orders of magnitude. 

5.2 Neutralino Induced Muon Fluxes 

It is often stated that one of the best indirect neutralino search strategies is to look for 
neutrinos produced by neutralino annihilations in the center of the Sun or of the Earth. 
Neutralinos from the galactic halo may get trapped into these astrophysical bodies, and 
begin to sink into the center, where the density enhancement produces annihilations into 
SM particles, among which neutrinos, the only particles which may later be detected. 
Neutrino telescopes with Km^ size may be able to single out the neutrino flux from 
the center of those celestial bodies, and distinguish it from the unavoidable neutrino 
background from cosmic rays interactions with the atmosphere. A crucial quantity in 
this game is the balance between the capture rate and the annihilation rate: if this is an 
equilibrium process, then the signal is at his maximum; otherwise large suppressions may 
occur. The equilibrium time scale inside the Sun is in most cases much smaller compared 
to that in the Earth, and though, depending on the SUSY model, equilibrium may not 
be reached even in the Sun, it very rarely occurs in the case of the Earth. As a result, 
the flux of muon neutrinos from the Earth is, in most cases, far below current and future 
sensitivities; for this reason, we do not show here results for this detection channel. 

We show in Fig. that the flux from the Sun may be large enough to be detectable 
at future experiments (we will use here the future sensitivity prospects for the IceCube 
experiment [02] )• A few models are even already excluded by current SuperKamiokande 
data |63]. It goes without saying that models with a larger annihilation cross section give 
larger rates: this point is clarified in the right panel of Fig. where we include also low 
relic density models. Interestingly, we find that, in any case, DM searches at neutrino 
telescopes will not probe neutralino masses larger than about 750 GeV. Notice that the 
overall gross features for spin- dependent rates and for the muon flux from the Sun are 
rather similar to each other: this does not come as a surprise, since the capture rate into 

®The same applies also for low relic density models, which we do not include in Fig. El 




Figure 14: Scatter plots of the expected muon flux at neutrino telescopes from neutralino annihi- 
lations in the center of the Sun, above a 1 GeV threshold. Models in (a) have relic abundances in 
the WMAP range, and are grouped according to the bino, higgsino and wino content of the light- 
est neutralino. In (b) we show instead a selection of all models grouped by relic abundance. We 
also show, respectively with a black and a green line, current and future experimental sensitivity 
limits. 
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Figure 15: A statistical analysis of the Visibility Ratio (signal-to-sensitivity) for the Muon flux 
from the Sun at Neutrino Telescopes for all models with a sufficiently low relic abundance (plots 
to the left) and for models with a relic abundance in the WMAP range (plots to the right). The 
two upper plots refer to the current exclusion limits, as in Fig. \14\ (a), black lines, while the two 
lower plots to the future exclusion limits. Fig. \14[ green lines. 



the Sun, mainly composed of nuclei with spin different from zero, depends in fact on the 
spin-dependent neutralino-nucleon cross section. 

The suitable neutralino candidate whose relic abundance lies within the WMAP range, 
and which will be detectable at neutrino telescopes, is a composite bino-higgsino state 
with a large enough spin-dependent cross section, and a mass below half a TeV. Models 
at larger masses (as pure Higgsinos and Winos with relic abundances within the WMAP 
range) will not be probed at future experiments. 

The statistical summary of our scan is shown in the histograms of Fig. ^1 Once 
again, current experiments only probe a marginal fraction of models, particularly if lying 
into the WMAP preferred range (histograms to the right). On the other hand, future 
experiments will be able to probe from 3 to 5% of the viable models. 

5.3 Correlating Direct Searches and Muon Fluxes 

In this section we analyze the correlation between neutralino direct searches and indirect 
searches at neutrino telescopes. In particular, as outlined above, we correlate, in Fig. ITBt 
spin-dependent neutralino-proton rates with the muon flux from the Sun, and the scalar 
neutralino-proton cross section with the rate of muons from the Earth (mainly composed 
by spin-less nuclei). As guidelines, we also include current and future sensitivities, at a 
putative neutralino mass of 1 TeV: models above (or at the right) the horizontal (respec- 
tively vertical) lines will be, or currently are, above projected, or current, sensitivity. The 
same color code as in Fig. has been used, although we plot here also models with low 
relic abundances. 




Figure 16: The correlation between the muon flux from the Earth and the Spin Independent 
neutralino-proton scattering cross section (a) and the correlation between the muon flux from 
the Sun and the Spin Dependent neutralino-proton scattering cross section (h). We show here 
a reduced sample of all models considered (including models with relic density below the WMAP 
range). Black circles represent models in which the lightest neutralino is predominantly a Bino, 
red squares represent Higgsinos and Blue diamonds Winos. As guidelines, we also include green 
and black lines representing future and current sensitivities at a neutralino mass of 1000 GeV. 



As a first remark, we point out the well known complementarity between direct and 
indirect searches: models which are not visible at spin-dependent searches will be ac- 
cessible at indirect searches in the muon-from-the-Sun channel, and, vice-versa, models 
with a too low muon flux from the Earth may well be above visibility threshold at direct 
spin-independent searches. 

The correlation between indirect searches at neutrino telescopes and direct detection 
experiments is clearly visible, though some comments are in order. First, binos always 
tend to have a smaller muon flux, at the same neutralino-nucleon cross section, than 
higgsinos or winos. This fact is due to the pair annihilation rate, which in the case 
of binos is typically more suppressed than for winos or higgsinos; the sufficiently low 
relic abundance of binos has been shown to occur in most cases thanks to coannihilation 
processes, which are however not present for pair annihilations in the potential wells of 
the Sun or of the Earth. Secondly, the correlation between (p^^^^^ and cx^p, though clearly 
present, is scattered over at least four orders of magnitude in (j)^^^^^ at a given a^p, for 
each neutralino type. Various factors contribute to this spread: the mentioned effect due 
to the pair annihilation rate, the details of the annihilation-capture interplay and possible 
enhancements due to kinematical effects for particular values of the neutralino mass. 

We recall that, if equilibrium is reached, both (p^^'^^^ and cr^p scale as the squared of 
the nucleon matrix element of the effective Lagrangian for the scalar neutralino-nucleus 
interaction |(£sc)P, while if it is not then [Hll 

0^^'-*^Oc|(Ac)r(^annt;)o (19) 

(c"ann^)o being the neutralino annihilation times the relative velocity in the zero velocity 
limit. In the case of the 0^"° - a^p correlation, we notice that large fluxes, corresponding 
to cases where annihilation and capture are in equilibrium, tend to have, as expected, an 
extremely strong correlation, which is lost when the signal is weaker, once again because 
equilibrium is not reached, and the dependence on (o"ann^')o again enters into the game. 
In particular, this is the case for binos, where coannihilation effects with a large variety 
of partners can drastically affect the actual neutralino annihilation rate with respect to 
what expected from cosmological abundance arguments: a coannihilating bino can in fact 
produce a sufficiently reduced relic abundance though featuring a large ((Tann^)o- 

5.4 Antiprotons, Positrons and Gamma Ray Searches 

Antimatter searches (see e.g. Ref. [HI]) have been recently shown to be an appealing search 
strategy for models with large annihilation rates, and therefore low relic abundances in 
the standard cosmological scenario [HS]. Large sources of uncertainties arise however 
in this context, mainly due to the computation of the secondary antiproton flux (the 
background), to the antimatter propagation in the Galaxy and in the Solar System 
and to uncertainties in the dark matter halo distribution in the Milky Way |i56j . 

We will here mainly refer to the analysis carried out in Ref. [HS]. The primary source 
diffusion in the Galaxy is modeled with an effective two-dimensional diffusion setup in 
the steady state approximation; the solar modulation effects are taken into account with 
the analytical one parameter force- fleld approximation of Ref. |SZ|. The secondary anti- 
matter fluxes are computed with the Galprop package, and give an excellent flt to the 
currently available data, respectively with a reduced = 0.82 for antiprotons and 0.95 
for positrons. 



In Ref. inS] a given SUSY model was ruled out if the sum of the primary and secondary 
antimatter fluxes gave a statistically unacceptable . For future experiments, a new pa- 
rameter 1^ was proposed, basically defined as the integral over antimatter kinetic energies 
of the ratio of the squared of the signal over the background. This quantity could then be 
compared against quantities related to the particular experimental facility at given data 
acquisition time scales, providing an assessment of the visibility of the given SUSY model. 
Here, for computational ease, we will take one representative energy bin, both at current 
and future antimatter search experiments, and compare the primary antimatter flux with 
the (current or projected) experimental sensitivity. Since the background is consistent 
with available data, we will indicate as visible those models where the antimatter flux of 
supersymmetric origin 0"susy gives a contribution which could have been, or which will 
be detected, i.e. 

0asuSY > 2 oTe^p, (20) 

where cxexp indicates the current, or projected, sensitivity in the given energy bin. We 
compared, for a sample of models, the method we use here with the fully statistical accu- 
rate strategy of Ref. [HS] , and we verified that the exclusion limits found were consistent 
with each other. 

The spectral features of primary antiprotons indicate a peak at energies of a few GeV, 
depending on the neutralino mass and composition. To this extent, we took the = 1.95 
GeV bin from the BESS-98 data [HHI, which turns out also to be very well fitted by the 
computed secondary fiux. As regards future perspectives, we took a putative Ep = 20 
GeV bin at the PAMELA experiment, for which we simulated the experimental sensitivity 
and the data binning after three years of data taking, following Ref. |^ : our choice is in 
this case motivated by the fact that space based experiments will take data in the large 
kinetic energy range, where balloon-borne experiments could not be sensitive enough, a 
region where the signal-to-background ratio is particularly large. We limited ourselves, 
however, to 20 GeV since our scan in principle allows for neutralino masses as low as 




Figure 17: Current and future discrimination sensitivities for antiproton fluxes, for a sample 
of models grouped by relic abundance. 
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Figure 18: A statistical analysis of the Visibility Ratio (primary flux over twice the error bars) 
for antiprotons with the Adiabatically contracted halo profile, for all models with a sufficiently 
low relic abundance (plots to the left) and for models with a relic abundance in the WMAP range 
(plots to the right). The two upper plots refer to the current data at Ep = 1.95 GeV, while the 
two lower plots to the future PAMELA projected sensitivity at Ep = 20 GeV. 
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Figure 19: The same statistical analysis for antiprotons, as in Fig. \18l but for the Burkert 
profile. 



around twice that value. For positrons, we chose the i?e+ = 11.09 GeV bin from the data 
of the HEAT-94/95 flights jTUI, and took = 23 GeV for the PAMELA experiment. 

Fig. El collects our results concerning antiprotons, in terms of the (j)^^'"^^ / {2aexp) 
ratio as a function of the neutralino mass, for current data (left) and at the forthcoming 
PAMELA experiment (right). We show here the results for the adiabatic halo profile. 
We also indicate the 5 — cr line, corresponding to a SUSY signal exceeding five times the 
current experimental sensitivity. The color code indicates different neutralino thermal 
relic abundances as computed in the standard cosmological scenario; as pointed out for 
three particular benchmark scenarios in Ref. [65j, only models with large annihilation rates 
can give measurable antimatter signals. In fact, models with thermal relic abundances 
n^/i^ < 0.01 are largely excluded by current data, and they will be thoroughly probed at 
future experiments. On the other hand, models within the WMAP range are not already 
constrained; only few of these models will be accessible in the future for neutralino masses 
up to 600-700 GeV. Notice the clear hierarchy of wino (upper line clustering) and higgsino- 
like (lower line clustering) neutralinos, which, due to the large branching ratio into gauge 
bosons final states, tend to have homogeneous antiproton fluxes, which only scale with the 
neutralino mass. This clearly allows model-independent predictions for these neutralino 
composition. Noticeably, a very large spread is instead present in the case of binos, 
where the dependence on the SUSY mass spectrum is critical for the computation of the 
antimatter flux. 

Our statistical analysis, for the antiproton fluxes, is shown in Fig. and Fig. 
respectively for the Burkert and for the adiabatically contracted proflle. With both halo 
proflles, models with large annihilation rates will be largely accessible to future space 
based experiments. On the other hand, models with thermal relic abundances within the 
WMAP range are not currently constrained, and will be one to two orders of magnitude 
below future sensitivity. We therefore stress that supersymmetric models yielding the 
correct dark matter thermal relic abundance in the standard cosmological scenario do not 




Figure 20: Current and future discrimination sensitivities for positron primary fluxes, for a 
sample of models grouped by relic abundance. 
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Figure 21: A statistical analysis of the Visibility Ratio (primary flux over twice the error bars) 
for positrons with the Adiabatically contracted halo profile, for all models with a sufficiently low 
relic abundance (plots to the left) and for models with a relic abundance in the WMAP range 
(plots to the right). The two upper plots refer to the current data at = 11.09 GeV, while 
the two lower plots to the future PAMELA projected sensitivity at = 23 GeV. 




Figure 22: The same statistical analysis /or positrons, as in Fia. \21V but for the Burkert profile. 



produce antimatter fluxes detectable at future space based experiments. 

The situation for positrons is rather similar to that of antiprotons, though, less promis- 
ing. The scatter plots of Fig. 1201 indicate that we will be able to detect, at best, models 
featuring a thermal relic abundance up to ^^^'^ ^ 0.01. The discrimination ratios for 
models within the WMAP range fall, in the best case scenario, approximately three or- 
ders of magnitude below future sensitivities. The histograms of Fig. 1^ and 1^ summarize 
the statistics we collected for the positron channel. 

To further detail on the correlation between thermal relic abundance and antimatter 
fluxes, and to compare antiproton and positron fluxes, we plot in Fig. |221 a sample of 
models, grouped by neutralino composition, in the 0"/(2cr)) plane. The case of 

antiprotons is showed in the left panel, while that of positrons in the right panel. First, 
notice that at a given Q^h"^, the discrimination ratio for antiprotons is everywhere one 
order of magnitude better than that of positrons. Secondly, a clear correlation between 
Q^h"^ and (j)°'/{2a) is present: Wino and higgsino like neutralinos lie on a relatively thin 
strip in this plane, which ends, in the WMAP preferred range, more than one order of 
magnitude below the visibility lines both for positrons and antiprotons detection. Binos 
scatter instead over a broad range of fluxes. In this respect, we generalize here the 
conclusions of Ref . , in that antimatter searches are a promising road towards indirect 
discovery of neutralinos provided some relic density enhancement mechanism is operative, 
and models with a large neutralino annihilation rate may account for the whole WMAP- 
inferred amount of dark matter. 

In the presence of a halo profile featuring a cusp in the Galactic Center (GC), a sizeable 
contribution to the 7-ray flux is expected to be provided by neutralino annihilations. 
The EGRET experiment, on the Compton Gamma- Ray Observatory, has reported a 
tentatively extended 7-ray source ~ 1.5° from the GC ^T], with a spectrum incompatible 
with the one expected from interactions of primary cosmic rays with the interstellar 




Figure 23: The correlation between the neutralino relic abundance and the projected sensitivities 
on antiproton (a) and positron (b) fluxes. Models are grouped according to the nature of the 
lightest neutralino. The WMAP range is indicated by a vertical green strip. 




Figure 24: A scatter plot for the ratio between the expected gamma ray flux from neutralino 
annihilations in the galactic center and the 2-a sensitivity in the bin at the largest energies 
sampled by EGRET. Models are grouped according to their relic abundance. 

medium, and compatible, on the other hand, with a WIMP-induced component )72j . 
Though ahernative explanations for this so-called EGRET excess have been put forward, 
these measurements may be turned around and used as a constraint, for a given halo 
profile, on the SUSY 7-ray contribution, and thus on SUSY models [HH]. The most 
stringent bound comes from the largest energy bin, with photon energies in the range 
4 GeV < -E^ < 10 GeV, where the standard background component is largely suppressed 
with respect to the measured 7-ray flux. In this respect, we directly compare the SUSY 
contribution with the 2 — a bound from the EGRET measurement, and therefore declare 
not viable (within a given halo profile) any SUSY model yielding a flux larger than the sum 
of the central value and twice the error bars. Our results, for the adiabatically contracted 
halo profile, are collected in Fig. |^ Once again we notice the usual ~ m^^ suppression, 
the clustering of the fluxes for wino and higgsino like neutralinos, and the correlation 
with the thermal relic abundance. Models within the WMAP range are not constrained by 
current data, while low relic density models with Q^h'^ < 0.01 are excluded. The statistical 
analysis reported in Fig. |2H1 highlights the large dependence of 7-ray results on the halo 
profile: assuming a cored profile (two lower histograms) this indirect detection channel 
fails to give any constraint on the SUSY parameter space, even in the SUSY models with 
the largest annihilation rate. Moreover (see the two histograms to the right) models in 
the WMAP range are quite far from being detectable, even in the most optimistic halo 
profile setup. Future measurements of the mentioned 7-ray source may either put more 
stringent constraint, or confirm alternative explanations for the EGRET excess |73j . 
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Figure 25: A statistical analysis of the Visibility Ratio (flux over twice the error bars) for 
gamma rays, for all models with a sufficiently low relic abundance (plots to the left) and for 
models with a relic abundance in the WMAP range (plots to the right). The two upper plots 
refer to the Adiabatically contracted profile, while the two lower plots to the Burkert profile. 

5.5 An Overview of SUSY dark matter Search Strategies 

In view of the previous considerations, we quantitatively summarize in Tab.|3]our results as 
regards the comparison of different supersymmetric dark matter detection strategies. We 
made use of two different, benchmark halo models, a cuspy profile {Adiabatically contracted 
profile) and a cored profile {Burkert profile), as outlined in jHEl- Results concerning 
direct detection and neutrino telescopes are very mildly affected by the halo model under 
consideration, and we therefore reported our results only for the Adiabatically contracted 
profile. 

In the column All Models we include both models whose thermal relic abundance falls 
in the WMAP preferred range and those with lower relic densities: we recall that we do 
not perform any neutralino density rescaling, under the hypothesis that some relic density 
enhancement or non-thermal production mechanism allows models with large annihilation 
rates to be compatible with the correct required amount of dark matter. 

Let us stress, first of all, that, in a statistical sense, the bulk of the supersymmetric 
parameter space compatible with the lower and upper WMAP bounds on the neutralino 
relic abundance has still not been probed by dark matter searches, as highlighted by the 
zeroes appearing in the third column. On the other hand, taking into account low relic 
density models, only antimatter searches are currently providing significant constraints 
on SUSY models, though with a large dependence on the assumed halo profile. On the 
other hand, neutrino telescopes and direct detection rule out only marginal portions of 
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Table 4: A summary of the statistical analysis of the various dark matter detection methods. 
The first two columns refer to all models (WMAP+low relic density), while the last two columns 
to models with relic abundance within the WMAP range only. For antimatter and gamma rays 
we indicate the results for both the Adiabatically contracted profile and for the Burkert profile. 

parameter space^. 

As regards future prospects, we find that a significant portion of the viable SUSY 
parameter space will be probed at large direct detection facilities (approximately 20%, 
quite independently of the relic abundance of the models). Perspectives at Neutrino 
Telescopes are less exciting, but we find that still around 5% of the parameter space will 
be accessible at IccCube. Antimatter searches will fail to provide any strong constraint 
on models with thermal relic abundance in the WMAP range; however, provided some 
relic density enhancement mechanism is operative, they could become a prominent road 
to SUSY dark matter discovery on future space based experiments (AMS, PAMELA). 
Antiprotons are found to be in any case a more promising dark matter indirect detection 
channel than positrons. 

6 Conclusions 

The present paper is devoted to a statistical analysis of the general minimal, flavor di- 
agonal, supersymmetric extension of the standard model. Resorting to a random scan of 
the MSSM i^arameter s])a.ee. in tlic soft breaking mass range between 50 GeV and 5 TeV, 

^This does not mean, of course, that there are no parameter space choices whose detection rates at 
direct searches and neutrino telescopes fall above current exclusion limits: what we find is that these 
choices are statistically marginal. 



we first studied the relative effectiveness of tlie various plienomenological constraints. We 
singled out 10^ viable models, also consistent with the upper WMAP bound, and per- 
formed a statistical study out of them, starting from the expected values of 6 ^ sj, rrih 
and {g — 2)^. We found that the CP even lightest Higgs boson mass is limited from above 
at about 132 GeV. The SUSY contributions to the muon anomalous magnetic moment 
have been shown to be always noticeably suppressed. 

We then turned to the topic of supersymmetric dark matter. We showed that the 
relic abundance of the viable models tend to be around one order of magnitude below the 
WMAP preferred range. This is due to the large number of models featuring a wino or a 
higgsino like lightest neutralino. If one requires the thermal relic abundance to be within 
the WMAP 2 — cr range, one finds that pure higgsino states masses cluster around 1 TeV, 
while winos around 1.6 TeV, killing any hope of detecting supersymmetric signals at the 
LHC. On the other hand, bino like neutralinos with the correct relic abundance feature, 
in the large majority of coannihilating partner, or a resonant annihilation cross 

section, and the resulting LSP mass spreads over a wide range. We analyzed the nature 
of coannihilating partners and the relative efficiency in the reduction of the bino relic 
abundance, and provided details on resonant bino annihilations into heavy Higgs bosons. 

Is it therefore correct to consider minimal supcrgravity as a benchmark MSSM sce- 
nario? The answer is probably not so straightforward. We found that models whose relic 
abundance lies in the WMAP range have in most instances a higgsino like neutralino, but 
that many models also feature a bino or a wino LSP: in this respect, mSUGRA introduces 
some kind of bias in predicting, over most of its parameter space, a bino-likc neutralino. 
As in mSUGRA, binos in the MSSM require the existence of relic density suppression 
mechanisms, like coannihilations or resonances: however, the general MSSM allows for 
a much richer plethora of coannihilating partner, and, consequently, a wider neutralino 
mass range. 

The second part of our analysis has been devoted to a comparison of supersymmetric 
detection strategies. We made use of two benchmark halo models, in which the density 
and velocity distributions have been self-consistently computed. The analysis has been 
carried out with the use of Visibility Ratios, which allowed an assessment of current and 
future detection channels and a comparison of different search strategies. 

We found that current data on dark matter searches provide weak constraints on the 
bulk of the general MSSM viable models. Future perspectives, for models whose thermal 
relic abundance lies within the WMAP range, appear to be quite promising for direct, spin- 
independent searches (we find that ~20% of viable models will be within future sensitivity) 
and, to a less extent, for neutrino telescopes looking at neutrinos produce by neutralino 
annihilations in the Sun. Antimatter searches were found to be largely dependent on the 
neutralino relic abundance: models within the WMAP range do not typically produce 
antimatter fluxes which will be detectable at future space based experiments. Taking into 
account the possibility of relic density enhancement mechanisms in the Early Universe, 
models with large annihilation rates could account for the whole of dark matter: in this 
case, antimatter searches put severe constraints on most SUSY models, and may become 
soon a prominent road towards dark matter discovery. This conclusion holds true quite 
independently of the details of the dark matter distribution in the inner part of our Galaxy, 
contrary to the case of gamma rays searches. Finally, we provided a summary table in 
which we quantify the statistical relevance of a given dark matter search strategy, both 
for low relic density models and for models with WMAP compatible neutralino thermal 



abundance. 
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